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Abstract 

We present a summary of a recent workshop held at Duke University on Partonic Transverse 
Momentum in Hadrons: Quark Spin-Orbit Correlations and Quark-Gluon Interactions. The trans- 
verse momentum dependent parton distribution functions (TMDs), parton-to-hadron fragmenta- 
tion functions, and multi-parton correlation functions, were discussed extensively at the Duke 
workshop. In this paper, we summarize first the theoretical issues concerning the study of partonic 
structure of hadrons at a future electron-ion collider (EIC) with emphasis on the TMDs. We then 
present simulation results on experimental studies of TMDs through measurements of single spin 
asymmetries (SSA) from semi-inclusive deep-inelastic scattering (SIDIS) processes with an EIC, 
and discuss the requirement of the detector for SIDIS measurements. The dynamics of parton 
correlations in the nucleon is further explored via a study of SSA in D (D) production at large 
transverse momenta with the aim of accessing the unexplored tri-gluon correlation functions. The 
workshop participants identified the SSA measurements in SIDIS as a golden program to study 
TMDs in both the sea and valence quark regions and to study the role of gluons, with the Sivers 
asymmetry measurements as examples. Such measurements will lead to major advancement in 
our understanding of TMDs in the valence quark region, and more importantly also allow for the 
investigation of TMDs in the unexplored sea quark region along with a study of their evolution. 

PACS numbers: 13.60.-r, 13.88. +c, 13.85.Ni 
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I. INTRODUCTION 



Understanding the internal structure of nucleon and nucleus in terms of quarks and 
gluons, the fundamental degrees of freedom of Quantum Chromodynamics (QCD) has been 
and still is the frontier of subatomic physics research. QCD as a theory of the strong 
interaction has been well-tested by observables with a large momentum transfer in high 
energy experiments. Our knowledge on the universal parton distribution functions (PDFs) 
and fragmentation functions (FFs), which connect the partonic dynamics to the observed 
hadrons, has been dramatically improved in recent years []] . As a probability density to find 
a parton (quark or gluon) inside a hadron with the parton carrying the hadron's longitudinal 
momentum fraction x, the PDFs have provided us with the nontrivial and quantitative 
information about the partonic structure of a hadron. 



In recent years, the hadronic physics community has extended its investigation of partonic 
structure of hadrons beyond the PDFs by exploring the parton's motion and its spatial 
distribution in the direction perpendicular to the parent hadron's momentum. Such effort is 



closely connected to the study and extraction of two new types o 
transverse momentum dependent parton distributions (TMDs 
parton distributions (GPDs) js 
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^parton distributions: the 
2|-l9|, and the generalized 
The ultimate knowledge of finding a single parton 
inside a hadron - involving both momentum and space information - could be encoded in 
the phase-space distributions of quantum mechanics, such as the Wigner quasi-probability 
distribution W(k,b), whose integration over the parton spatial dependence (6) leads to the 
TMDs, while its integration over transverse momentum (k) provides the parton's spatial 
distribution that is relevant to the GPDs. A quantum field theory version of the phase- 
space distributions, in terms of the matrix element of the Wigner operator, was discussed 
in Ref . 17| • Understanding both the momentum and spatial distribution of a parton inside 



a hadron in terms of the more general Wigner distributions could be the central object of 
future studies on partonic structure. Knowledge of TMDs is also crucial for understanding 
some novel phenomena in high energy hadronic scattering processes, such as, the single 



transverse spin asymmetries 



18H24| and small- x saturation phenomena 25l-l32|. 



II. RECENT THEORETICAL DEVELOPMENT ON TMDS AND EXPERIMEN- 
TAL ACCESS 



Like the PDFs, the TMDs and GPDs carry rich information on hadron's partonic struc- 
ture, while they are not direct physical observables due to the color confinement of QCD 



dynamics. It is the leading power QCD collinear factorization theorem [33| that connects the 
PDFs to the hadronic cross sections with large momentum transfers: Q's S> Aqcd- in order 
to study the TMDs, we need the corresponding TMD factorization theorem for physical ob- 
servables that are sensitive to parton's transverse motion and the TMDs. Such observables 
often involve two very different momentum scales: Qi ^> Q 2 ^ Aq CD , where the large Qi 
is necessary to ensure any perturbative QCD calculation while the small scale Q2 is needed 
so that these observables are sensitive to the parton's transverse motion. The transverse 
momentum distribution of single hadron production in semi-inclusive deep inelastic lepton- 
hadron scattering (SIDIS) and Drell-Yan lepton pair production in hadronic collisions are 
two natural examples. The TMD factorization for these two processes have been carefully 
examined 
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361] . However, the TMD factorization in QCD is much more restrictive than 
the leading power collinear factorization. The conventional TMD factorization works for 
three-types of observables with only two identified hadrons: single hadron px distribution 
in SIDIS, the px distribution of Drell-Yan type process, and two-hadron momentum imbal- 
ance in e + e~ collisions. But it has been shown to fail for observables with more than two 



dentified hadrons 37H41] 



Important aspects of the TMD parton distributions, such as the gauge invariance, th e 

ro 1 
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e of gauge links, and the universality, have been explored in recent years |20H24 
44j. Like the PDFs, the definition of TMDs is closely connected to the factorization of 
physical cross sections, and it is necessary for the TMDs to include all leading power long- 
distance contributions to the physical cross sections if they could be factorized. All leading 
power collinear gluon interactions are summed into the gauge links in the definition of 
the TMDs. It is the gauge link that makes the TMDs gauge invariant and provides the 
necessary phase for generating a sizable single transverse spin asymmetry (SSA) in SIDIS 



and Drell-Yan processes 18l-l23|. However, unlike the PDFs, which are universal, the TMDs 
could be process dependent due to the fact that the initial-state and final-state collinear 
gluon interactions are summed into two different gauge links. That is, the TMDs extracted 



from SIDIS could be different from those extracted from Drell-Yan processes because of the 
difference in gauge links. Although the TMDs are not in general universal, it could be shown 
from the parity and time-reversal invariance of QCD dynamics that the process dependence 



of the spin-averaged as well as spin-dependent TMDs is on 
as the parity and time-reversal modified universality 
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y a sign, which was referred to 



24] . An important example of the 



modified universality is that the Sivers function extracted from the SIDIS measurements is 
opposite in sign from the Sivers function extracted from the Drell-Yan process. The test of 
the sign change of the Sivers function from SIDIS to Drell-Yan is a critical test of the TMD 
factorization. 

r~| 

At leading twist there are eight TMD quark distributions [9]: three of them, the unpo- 
larized, the helicity and the transversity distributions, survive in the collinear limit, while 
the other five vanish in such a limit. Each TMD quark distribution explores one unique 
feature of the quark inside a polarized or unpolarized nucleon. For example, the Sivers 
function 3), l^j] provides the number density of unpolarized partons inside a transversely po- 
larized proton, while the Boer- Mulders function |6| gives the number density of transversely 
polarized quarks inside an unpolarized proton. Although we have gained a lot of informa- 
tion on the collinear PDFs and helicity distributions, we know very little about quark's and 
gluon's intrinsic transverse motion inside a nucleon. Recent measurements of multiplicities 
and double spin asymmetries as a function of the final transverse momentum of pions in 



SIDIS at JLab 



46 



47] suggest that transverse momentum distributions may depend on the 



polarization of quarks and possibly also on their flavor. Calculations of transverse momen- 
tum dependence of TMDs in different models 
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5 11 ] and on lattice [52j, |53[ indicate that 



dependence of transverse momentum distributions on the quark polarization and flavor may 
be very significant. 

Among the TMDs vanishing in the collinear limit, the Sivers function is the best known 



and has been phenomenologically extracted by several 



azimuthal distribution of a single hadron in SIDIS 



groups mainly from analyzing the 



541457] . However, in the case of posi- 



tive hadrons, where a signal has been seen, the measurements of HERMES [58] and COM- 



PASS 



59J experiments are only marginally compatible: the asymmetries measured by COM- 



PASS are somewhat smaller, and seem to indicate an unexpected dependence on W, the mass 
of the hadronic final state. For the transversity distribution, there is only one phenomeno- 
logical extraction by combining the SIDIS and the e + e~ data 6CH62]. and information on 



the rest of the TMDs is rather scarce. Nevertheless, these recent results have already gen- 
erated great excitement, which is evident from the increasingly active theoretical activities, 
including modeling and lattice QCD calculations, and planning of future experiments. 



A number of experimental facilities, such as COMPASS 63[ at CERN, CEBAF with its 
12 GeV upgrade at Jefferson Lab, RHIC at Brookhaven National Lab, Belle at KEK, and in 
particular, the planned Electron-Ion Collider (EIC), will play a complementary but crucial 
role in determining these TMD parton distributions. Among three types of processes where 
the TMD factorization could be valid, the SIDIS might be the best place to study the TMD 
parton distributions because of the easy separation of various TMDs with well-determined 
distributions in the azimuthal angles between the spin, the leptonic plane and the hadronic 
plane, in addition to the much higher event rates. With a broad energy range and a high 
luminosity, the future EIC will be an ideal place to extract the TMDs in a multi-dimensional 
phase space with a high precision. Precise measurements of these new distributions could 
provide us much needed information on the partonic structure of nucleon (nucleus) in order 
to address the fundamental questions concerning the decomposition of the nucleon spin, and 
the QCD dynamics responsible for the structure of the nucleon. 

The transverse momentum dependence could also be introduced to the hadronization 
process to get the TMD parton-to-hadron fragmentation functions. For a quark to fragment 
into a spinless hadron, such as a pion, there are only two possible fragmentation functions 
at the leading twist: the unpolarized fragmentation function and the Collins function 64|, 



which is responsible for generating azimuthal asymmetric distribution of hadrons from the 
hadronization of a transversely polarized quark. The Collins function has been extracted 
from recent experiments (Belle [65], HERMES 66|, |67j, and COMPASS 68]), and was found 
to be nonzero. Precise measurements of TMD fragmentation functions provide a new window 
to explore the dynamics of hadronization. 

For cross sections with one large momentum transfer, or several momentum transfers at 
the same scale, it is more natural and appropriate to use the collinear factorization approach 
and to expand the cross sections as an inverse power of the large momentum transfer. Al- 
though the leading power term dominates the contribution to the cross sections, it does not 
contribute to the SSA, which is proportional to the difference of two cross sections with the 
spin vector reversed. Like the TMD factorization approach, the SSA in the collinear factor- 
ization approach is also generated by the active parton's transverse motion, but, as a net 



effect after integrating over all possible values of transverse momentum. Within the collinear 
factorization approach, the SSA is effectively generated by the quantum interference of two 
scattering amplitudes: a real amplitude with one active parton and an imaginary part of 



an amplitude with an active two-parton composite state [69h74|. The QCD quantum inter- 
ference between the amplitude to find a single active parton and that to find a two-parton 
composite state is represented by a set of new twist-3 three-parton correlation functions. 
Unlike the PDFs, which have the probabilistic interpretation of number densities to find 
a parton within a hadron, these new three-parton correlation functions provides the direct 
information on the strength of color Lorentz force and/or magnetic force inside a spinning 
proton. The twist-3 contributions are accessible in various spin-azimuthal asymmetries in 
SIDIS depending on the helicity of the lepton or the hadron. Significant higher-twist asym- 
metries have been reported by the HERMES 75H77| and COMPASS Collaborations 78j as 



well as the CLAS and Hall-C Collaborations at JLab [46|, |47j, 179[ . Higher- twist observables, 
such as longitudinally polarized beam or target SSAs, are important for understanding long- 
range quark-gluon dynamics, and the future EIC due to the wide range in Q 2 , will be an 
ideal place to pin them down. 

Both the TMD factorization approach and the collinear factorization approach at twist-3 
provide a viable mechanism to generate the SSA, but, with a very different physical picture. 
This is because they cover the SSA in two very different kinematic regimes: Q\ ^> Q2 ^ 
Aqcd for the TMD approach while Qi ^> Aqcd with i = 1,2, ... for the twist-3 approach. 
Further study has shown that the TMD approach is consistent with the twist-3 approach for 
ihe SSA phenomena in a perturbative region, Qi ^> Q 2 Aqcd, where they are both valid 
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More recently, the evolution equations for the transverse momentum moments 
of these TMDs have also been investigated, which opens a path for the systematic QCD 
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calculations of SSA beyond the leading order in a s 

Like TMD quark distributions, we could also construct TMD gluon distributions. But, 
unlike the quark, gluon does not interact with any colorless particles at the lowest order. 
Due to the restriction on the color flow for the TMD factorization to be valid, we only 
have very few observables that might give the direct access to the TMD gluon distributions, 
such as the Higgs production at low pt with an effective gg — > H° vertex, the momentum 
imbalance of two isolated photons via an effective gg — > 77 vertex, and back-to-back jets or 
heavy quark pair production 3, 0-911 in ep and in pp or pp collisions. 



On the other hand, many more observables could access the gluonic sector of twist-3 
approach to the SSA. Heavy flavor production in the DIS regime is a direct probe of gluon 
content of the colliding hadron or the nucleus. In particular, the SSAs of open-flavor (anti)D 
(or B) meson production in the DIS regime provides a unique opportunity to measure tri- 
gluon correlation functions [92] , which are closely connected to the gluon's transverse motion 
and color coherence inside a transversely polarized nucleon. A more recent study shows that 
there are four tri-gluon correlation functions |93|. The co-existence of these different tri- 
gluon correlation functions, which represent the long-distance quantum interference between 
a gluon state and a two-gluon composite state, is a unique feature of the non-Abelian color 
interaction. Motivated by recent calculations [92] , preliminary simulations of the event rate 
and asymmetries for some realistic EIC energies and possible detector coverage have been 
carried out, which will be presented in the later section of this paper. 

Like the PDFs, the TMDs and the GPDs are non-perturbative functions and should be 
extracted from the experimental measurements of cross sections or asymmetries in terms of 
relevant factorization formalisms. In order to get a better picture of the proton's partonic 
structure from the limited information extracted from the PDFs, the TMDs, and the GPDs, 
model calculations of these distributions are valuable. There have been many interesting 
model studies recently, see for example 50, 90, 94-9^]. These models and their calculations 



could play a very important role as a first step to describe the experimental observations, 
to give an intuitive way to connect the physical observables to the partonic dynamics, and 
to provide key inputs to the partonic structure of the nucleon, which will help us to address 
the fundamental questions, such as how the quark spin and its orbital angular momentum 
contribute to the nucleon spin? 

More importantly, very exciting results of TMDs have come from the Lattice QCD calcula- 



tions recently 52, 



53 



98( , indicating that spin-orbit correlations could change the transverse 
momentum distributions of partons. Notable results from Lattice QCD have been obtained 
for the impact parameter dependent parton distributions, which have a close relation to some 



interesting TMDs [99j. With the improvement of computer speed and simulation algorithms, 
more and more accurate results on the partonic structure from Lattice QCD calculations 
will become available soon. 

The experimental investigation of multi-dimensional spatial distributions of a parton (or 
color) inside a bound proton, in terms of the TMDs and the GPDs or the "mother" Wigner 
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distributions has just started recently. Future machines, like EIC, could supply high quality 
data by scattering polarized leptons off polarized nucleons. For semi-inclusive reactions, the 
data with large Q 2 and small Pt are dominated by the TMDs, while data with large Q 2 and 
large Pt have the most contributions from pQCD corrections convoluted with collinear PDFs 
or multiparton correlation functions. The transition from one regime to the other might be 
the most interesting aspect and should be carefully studied. The investigated x-region should 
be as wide as possible to cover both the valence quark region and the unexplored sea quark 
region. 

In summary, while there has been progress on several fronts in the theoretical devel- 
opments for understanding the transverse momentum dependent parton distributions and 
fragmentation functions, it is just a beginning for us to explore the full picture of partonic 
structure inside a nucleon. The future Electron-Ion Collider is a much needed machine to 
probe the partonic structure of a bound nucleon, to quantify the role of gluons and the color, 
and to help approach the fundamental question of strong interaction - the confinement of 
the color. Given below is a list of questions which have been discussed in various meetings 
in connection with the TMDs: 



Q 2 evolution. The transverse momentum dependence of the TMDs certainly depends 
on the large scale Q where the TMDs were probed. Al though th e energy evolution 



equation has been derived for the TMD distr ibuti ons 



I ()Q|- 102j |. very few explicit 



calculations have been performed to date (e.g. [1Q3() to study the Q 2 dependence of 
the associated experimental observables, such as the azimuthal asymmetries. This Q 2 
evolution is not only an important theoretical question, but also a crucial point to 
investigate experimentally. An EIC machine with a wide range of coverage on x and 
Q 2 for SIDIS processes will provide a great opportunity to study the scale dependencies 
of the TMDs in detail. 

The transverse momentum distribution of the TMD observables in principle, has three 
characteristic regions: intrinsic, resummation, and perturbative. In practice, a Gaus- 
sian transverse momentum distribution has been used to fit the existing experimental 
data in order to extract the TMDs. At low collision energy, there is not much phase 
space for the gluon shower around the hard collision, the active parton's intrinsic pt 
distribution dominates, and therefore, a Gaussian distribution should be a good ap- 



proximation. However, with a much higher collision energy at an EIC machine and 
a larger phase space for the gluon shower, the Gaussian distribution will not be ade- 
quate to describe the observed distribution in SIDIS. Instead, a distribution with a 
proper r esummat ion of large logarithmic contributions from the gluon shower should 
be used 
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1061 ] . When is as la rge as the hard scale Q, a perturbative calculated 



Pt should be more relevant 
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106| . Investigation of the resummation or "matching" 



region, especially as a function of Q 2 , will provide an important test of the theoretical 
framework, i.e. TMD factorization. 

Relation of TMDs to the parton orbital angular momentum. There have been qual- 
itative suggestions abou t the connection of the quark orbital angular momentum to 
the TMDs 
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107] . However, we still do not have a rigorous way to build this 



connection. Certainly, model calculations will help shed light on this important issue. 

Besides the connection to the parton's orbital angular momentum, the TMDs should 
provide much richer information on nucleon structure in momentum space. We need 
more theoretical investigations along this direction. 

Global study at next-to-leading order. So far, all phenomenological studies are limited 
to the leading order in perturbative QCD. We have to go beyond this simple picture 
to build a systematic framework to extract the TMDs. 

Small-x parton distributions. The investigation of TMDs at small- x (sea) has started, 



and progresses have been made recently 



29 



32 



108 



109| . In particular, it was found 



that di-hadron/dijet correlation in DIS processes in eA collisions can be used to probe 
the Weiza cker- Willi ams gluon distribution formulated in the color-glass-condensate 
formalism 



108 



109J. However, more theoretical studies are needed to build a rigorous 



connection between the TMD and CGC approaches. 

Universality of the TMDs. Much of the predictive power of QCD factorization relies 
on the universality of nonperturbative distributions. More work is needed to better 
understand the process dependence of the TMDs and their connections to what Lattice 
QCD can calculate, which is crucial for the predictive power of the TMD factorization 
and physical interpretation of the TMDs. 
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We believe that the TMD community as it addresses the above and other important 
questions will naturally make the case for an EIC stronger and will be ready for the new 
era of QCD and hadron structure. In the next section, we will present simulations that 
have been carried out with a goal of addressing the aforementioned questions with a high 
luminosity electron-ion collider. 



III. SIDIS AT AN EIC 



A. Kinematics 



In an EIC, a beam of electrons collides with a beam of ions. The SIDIS process requires 
to detect both the scattered electron and one of the leading hadrons produced in the final 
state. In general, the process can be expressed as: 

1{PI) + N(P) -> £'(P e ) + h(P h ) + X (1) 

where £, N, £' and h denote the initial electron, the initial proton 0, the scattered electron, 
and the produced hadron in the final state, respectively. All the four-momenta are given in 
parentheses. 

Under the one-photon exchange approximation, the four-momentum of the virtual photon 
is expressed as q = P % e —P e and the four momentum transfer square is q 2 = —Q 2 . The relevant 
Lorentz invariant variables are defined as: 

x = Sr q >y = ^ z = ^ s = & + P?*^-Ep. (2) 

Here, x, also referred to as Bjorken x, represents the initial nucleon momentum fraction 
carried by the parton in the infinite momentum frame, y and z are the fractional momentum 
carried by the virtual photon and the leading hadron, respectively, and s is the center-of- 
mass energy squared of the initial electron-nucleon system. The last approximation in Eq. [2] 
is made by neglecting the masses of the electron and the nucleon, which are much smaller 
than the center-of-mass energy at EIC kinematics. 
With approximations, one can immediately obtain 

Q 2 = x -y- s, (3) 



1 We assume that the ion is proton for simplicity. In principle, the same argument also applies to any ion 
beam. 



11 



FIG. 1. Definitions of azimuthal angles 4>h and (frs, and the hadron transverse momentum for SIDIS 
in the ion-at-rest frame j^]. 

which clearly illustrates the relation between x and Q 2 at fixed s. 

In addition to the aforementioned Lorentz invariant variables, there are a few frame- 
dependent kinematic variables, 4>s, 4>h, and Pt (the target spin angle, the azimuthal angle 
and the transverse momentum of the leading hadron), which are also essential to SIDIS 
process. They are defined according to the Trento convention as illustrated in Fig. [T] in the 
nucleon-at-rest frame 0. 

B. Phase Space Coverage 



In this section, we discuss the SIDIS phase space coverage mainly with the 11+60 GeV 
configuration, which represents all GeV electron beam colliding with a 60 GeV proton 
beam. In the simulation, the scattered electrons are generated in momentum P e > 0.7 
GeV/c, polar angle 2.5° < 9 e < 150° and full azimuthal angle. Fig. [2] shows the Q 2 
vs. x phase space for the 11+60 GeV configuration. Since we are mainly interested in 
the DIS region, the following cuts, Q 2 > 1 GeV 2 and W > 2.3 GeV, are applied. In 
addition, the 0.05 < y < 0.8 cut is also applied. Here, the 0.8 cut-off is chosen to reflect the 
lowest detectable energy of the scattered electrons, which is usually limited by the hardware 

2 More generally, the <f>s, 4>h : and Pt are defined in the collinear frame, where the virtual photon moves 
col-linearly with the initial nucleon. The nucleon-at-rest frame is a special situation of the collinear frame. 
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FIG. 2. Phase space of Q 2 vs. x with DIS FIG. 3. Momentum vs. polar angle in the lab 
and y cuts illustrated. No SIDIS cross section frame for the scattered electron after weighting 
is applied. The "p" in the legend refers to the by the SIDIS differential cross sections. Here, 



acceptance and uncertainties in the radiative correction The 0.05 cut-off of y is limited 
by the resolution of x. As shown in the following equation: 



assuming a fixed momentum resolution SP e /P e , the resolution of x increases dramatically 
at small y. Therefore, the y > 0.05 cut is applied in the simulation in order to maintain a 
reasonable x resolution for forward electron detection. With the above cuts applied, Fig. [3] 
shows the distribution of momentum vs. polar angle of the scattered electron in the lab 
frame after weighting each event by the SIDIS differential cross section. Two observations 
are made to reflect the needs in the detector design: 

• Most of the scattered electron are concentrated in the high-momentum region (closer 
to the initial electron momentum), which is corresponding to the small y region. This 
is consistent with the fact that SIDIS cross section are larger at smaller y. 

3 The larger the y value is, the more radiative correction should be applied, which would lead to larger 
systematic uncertainties. 



fact that the ion beam is the proton beam. 



0° represents the momentum direction of the 



initial electron beam. 
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P T <1.0 GeV 




FIG. 4. Momenta vs. polar angles for the detected hadron in the lab frame (weighted by differential 
cross section). The 180° represents the initial momentum direction of the ion beam. 



• No electrons are distributed at very forward angles (< 5°) due to the Q 2 > 1 GeV 2 
cut applied. There is no essential need to detect very forward-angle electrons. 

For SIDIS process, more cuts are applied on the hadron side. They are 0.2 < z < 0.8 
and Mx > 1.6 GeV cut, where Mx is the missing mass of the X system in Eq. [TJ The 
0.8 cut-off in z excludes the events from exclusive channels. The 0.2 cut-off is required to 
stay in the current fragmentation region, where the detected hadron can be used to tag the 
struck quark. In addition, we also apply a low Pt cut (Pt < 1 GeV/c) for the TMD physics, 
and a Pj- > 1 GeV/c cut for the large Pt physics. Fig. H] shows the momenta of detected 
hadrons vs. polar angles in the lab frame. Events are weighted by the SIDIS differential 
cross section. In this simulation, the hadrons are generated for 0.7 GeV/c < Phadron < 
10 GeV/c, full polar and azimuthal angular coverages. The Pt < 1 GeV/c cut is applied. 
Three observations are made to reflect the needs in the detector design: 

• Most of the hadron events are concentrated in the momentum region of 0.7-7 GeV/c. 
There is no essential need to cover very high momentum region. 
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P T <1 .0 GeV | 




x 

FIG. 5. Mapping of SIDIS phase space of dif- 
ferent energy configurations with proton beam. 
The 12 GeV phase space is shown in the black 
band. The blue, red, green bands represent the 
phase space of 3+20, 11+60 and 11+100 GeV 
configurations, respectively. 



P T <1 .0 GeV | 




x 

FIG. 6. Mapping of SIDIS phase space of differ- 
ent ion beams, given the fixed accelerator con- 
figuration. 



The hadrons have a wide distribution of the polar angle in the lab frame. 



No essential need to cover the very backward angle for the hadron lab polar angle. 
However, a large backward angular cover age i s important for the study of SSA from 



SIDIS in the target fragmentation region 



The upcoming JLab 12 GeV upgrade wou 



110]. 



and high x region due to the smaller s 



111 



d a ccess the SIDIS phase space at low Q 7 



- 1117] . The blac k band i n Fig. shows the 



111 



1 121 ] . In order to bridge 



phase space of the approved 11-GeV SoLID SIDIS experiment 
between the phase spaces of 11+60 GeV configuration and the JLab 12 GeV upgrade, a 
low-energy configuration of EIC, e.g. 3+20 GeV configuration is strongly desired. Such a 
configuration would overlap with both phase spaces of the 11+60 GeV configuration and the 
JLab 11-GeV fixed-target experiment. In addition, a higher-energy configuration, 11+100 
GeV (green band in Fig. E]), would extend the study of SIDIS process to even lower x and 
higher Q 2 regions. 

In order to achieve a quark flavor separation from the SIDIS data, measurements with 
both proton and neutron are essential. Since there is no free high-energy, high-intensity 
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P T (GeV) (ion-at-rest frame) P T (GeV) (ion-at-rest frame) 



FIG. 7. The detected momenta of hadrons vs. Pt (left) and polar angles vs. Pt (right) at Pt >1 
GeV/c. 

neutron beam available, one has to use the light ion beam instead. A deuteron beam is a 
natural choice. In the polarized case, the 3 He ion has a unique advantage as the effective 
polarized neutron beam; the ground state of 3 He is dominated by the S-state, where the 
two protons are arranged with spin anti-parallel to each other. Therefore, the 3 He spin 
is dominated by the neutron spin. However, the phase space of the ion is not the same 
as that of the proton. Given a fixed accelerator configuration, momentum per nucleon 
in an ion is proportional to Z/A, in which Z is the atomic number, and A is the mass 
number. Therefore, the light ion beam would lead to a smaller s with the same accelerator 
configuration. Fig. [6] illustrates different mapping of these three ion beams (accelerator: 
11+60 GeV configuration 0) . Therefore, 

• The lowest achievable x value for quark flavor separation is limited by the light ion 
beam rather than the proton beam. 

• The highest achievable Q 2 value for quark flavor separation is also limited by the light 
ion beam. 

At high Pt region (Pt > 1 GeV/c), the requirements on the hadron detection are shown 
4 60 GeV represents the momentum for proton. 
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in Fig. [7J The momenta of the hadron (left) and the lab polar angles of the hadron (right) 
are plotted vs. Pt- We make the following observations: 

• The hadron momentum range will increase with the increment of Pt- 

• The hadron lab angles distribute widely over the entire phase space. 

• It is not essential to cover the very backward angular range for the hadron lab polar 
angle. 



C. Transverse Single Spin Asymmetry Measurements for Light Mesons 

At an EIC, the transverse single spin asymmetry (TSSA) measurements with an unpo- 
larized electron beam and a transversely polarized proton (or effective neutron) beam can 
provide rich information on the transverse spin structure of the nucleon. Three leading twist 
TMDs, transversity, Sivers and pretzelosity distributions can be accessed. A large Q 2 cover- 
age of EIC will allow for a detailed study of the Q 2 evolution of the TMDs. The coverage in 
the small x region is essential to study sea quark TMDs. In particular, the light-meson (tt ± , 
SIDIS process will allow for a map of the TSSA for the sea quarks at low Q 2 and for 
the valence quarks at high Q 2 . Since TSSAs of SIDIS depend on four kinematic variables 
x, z, Q 2 , and Pt, a complete understanding would require mapping the TSSA in 4-D phase 
space. Therefore, a high luminosity machine is essential. In this section, we will illustrate 
the impact of a high luminosity EIC on TSSA measurements. 



1. Monte-Carlo Method 

Since most of the TSSAs are relatively small @, the projected uncertainty of the measured 
asymmetries can be approximated as: 

1 1 



* A *= PPP f 'TT'^ 

1 1 (5) 

PePlP N f D VK a iv 



5 As shown later, the expected asymmetries on proton is about a few percent, while the asymmetries on 
light ion are even smaller due to dilutions from spectator nucleon(s). 
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where P e , Pi and Pn are the polarizations of the electron, ion beam, and effective polarization 
of the nucleon. The fo is the effective dilution factor, and N raw is the raw measured counts 
summing over the two spin states. In the case of a proton beam, P^ = 1 and fr> = 1. In 
the case of a 3 He beam, P N = 87.5% and fjj ~ 0.3. 

In addition, N raw is the measured counts summing over the two spin states. Therefore, 
it is only proportional to the unpolarized cross section. The following Monte-Carlo proce- 
dure is adopted to obtain the projections on the separated Collins, Sivers and pretzelosity 
asymmetries. 

• Simulate scattered electrons and generate pions/kaons uniformly in both momentum 
and coordinate space in the lab frame. Cuts are applied to mimic the expected detector 
acceptance. 

• Apply the SIDIS cuts as described in Sec. IIIIBI 

• Calculate the SIDIS differential cross section for each accepte d eve nt. In this step, 



one has to calculate a 6x6 Jacobian matrix (see Sec 9.7 of Ref. 



118 



for the complete 



derivation) to transform the SIDIS differential cross section to the lab frame. 

Combining with the expected luminosity, running time, one can calculate the expected 
raw number of events in each of the 4-D kinematic bin. 

The projected uncertainties on the raw asymmetry are obtained after including the 
beam polarizations, nucleon effective polarization and the effective dilution factor. 

Additional factors are introduced to mimic the increase of uncertainties due to the 
azimuthal angular separation of Collins, Sivers, and pretzelosity asym metri es. The 
detailed discussion of these factors can be found in Appendix II of Ref. 11 11 ] . In the 



case of a full and uniform azimuthal angular coverage of 4>h and <ps, these factors equal 
to \/2, and are independent of the number of terms used in the fitting. More generally, 
they depends not only on the angular coverage, but also on the event distribution of 
the azimuthal angle. The factor on the Collins asymmetry is the same as that of the 
pretzelosity asymmetry, and slightly different from that of the Sivers asymmetry. In 
practice, these factors are calculated based on the simulated event distribution from 
Monte-Carlo. 
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2. Calculation of SIDIS Differential Cross Section 



In this section, we provide more details in how the SIDIS differential cross sections are 
calculated. The cross section at small Pt (Pt <1 GeV/c) is calculated based on the TMD 
formalism JjJ. 

da a 2 y 2 / 7 2 \ f ^ 

dxdydipdzd<p h dP^ xyQ 2 2(l—e)\ 2x J C/f/,T ' 

where 

2Mx 

and a is the fine structure constant. The angle ip is the azimuthal angle of I' around the 
lepton beam axis with respect to an arbitrarily fixed direction, which in case of a transversely 
polarized ion, it is chos en to be the direction of S. The corresponding relation between ip 



and <ps is given in Ref. 119} ; in deep inelastic kinematics one has dip ~ dips- The structure 
function, Fuu,t on the r.h.s. can be expressed as: 

F UUtT = xJ2 e l [ d 2 T f(x,p 2 T )D(z,\P T -zp T \ 2 ), (8) 

which depends on x, Q 2 , z and P%. Here, the first and second subscript of the above structure 
function indicate the respective polarization of the lepton and the ion beam, whereas the 
third subscript "T" specifies the polarization of the virtual photon with respect to the virtual 
photon momentum direction. The conversion to the experimentally relevant longitudi nal o r 



transverse polarization w.r.t. the lepton beam direction is straightforward and given in 119j . 
The ratio e of longitudinal and transverse virtual photon flux in Eq. [6] is given by 

1„2 i U2„2' W 



l-y + ±y 2 + 



v 



In order to calculate Fuu,t> the Gaussian ansatz [9i. lll9j |. is adopted for the transverse mo- 
mentum dependent parton distribution f{x,p^) (TMD) and fragmentation function D(x, k 2 ) 
(FF): 

f(x,p 2 T ) = f(x,0)exp(-R 2 H p 2 T ) (10) 
D(z, k 2 T ) = D(z, 0) exp(-i^). (11) 



Therefore, Fjjut becomes 



E<:/(")"(^, (12) 



'UU,T = X 

z 

a 



19 



FIG. 8. Feynman diagrams corresponding to the NLO contribution at high Pt region. 



where G(Qt] R) = (-R 2 / 71 ") exp(—Q T R 2 ), i.e., a Gaussian of which Qt = Pt/z, and the fall- 
off is determined by a radius R. Such radius is related to the radii Rh and Rh governin g the 
fall-off of f(x,p 2 T ) and D(x, k 2 T ) as R 2 = R\ R 2 h /(R 2 H + R 2 h ). Furthermore, the CTEQ6M [120 1 
is used to parametrize the parton distribution function (PDF) f(x, 0). The parametriza- 
tion of the unpolarized fragmentation function D(z,0) is from Ref. 1211 ] . R\ and R\ are 
assumed to be 0.25 GeV 2 and 0.2 GeV 2 [54J], respectively. While the assumption of the x 



and p\ factorization in addition to the Gaussian ansatz for the p\ dependence has been used 
widely in the literature, a statistical model for TMDs 51[ has been developed recently which 



involves a non factorisable x and p 2 T dependence, and the comparison wit h those w hich also 



have non factorisable TMDs, based on the relativistic covariant method 
made. We remark that in no TMD model considered so far 



50 



122- 



20, 195H93, 



24 1 has been 



125 



133fl the 



Gaussian Wz is strictly supported, although so m e support it approxunately W and so 
does phenomenology [94J |. 

The cross section at larg e Pt (Pt >1 GeV/c) is dominated by the pQCD higher-order 



collinear contributions 



134j. The dominant partonic processes are shown by the Feynman 



diagrams in Fig. [HJ The quark can emit a hard gluon, or be generated by the gluon through 
pair production. The expression of the Q(a s ) (LO ) calculation of the differential cross section 



134j, in which the R 2 H and R\ are assumed to 



in the high region can be foun d in Ref. 
be 0.28 GeV 2 and 0.25 GeV 2 [l34], respectively. 

The cross sections in both low and high Pt regio ns are expressed as 



da 



dx dy dip dz d4>h dPii 

in the ion-at-rest frame. A 6x6 Jacobian matrix |118| is derived to convert the calculated 
differential cross section to the lab frame in terms of the kinematic variables of the final-state 
lepton and hadron: P e , 9 e , <f> e , P^, Oh, 4>h- 

We have carried o ut s tudies to compare results from our approach described above with 



those from EICDIS 



1351 ] . which is based on PEPSI generator, and results from PYTHIA 
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Ion 


11+60 GeV 


3+20 GeV 


11+100 GeV 


Polarization 


P 


9.3 x 10 40 cm" 2 


3.1 x 10 40 cm~ 2 


3.1 x 10 40 cm" 2 


1 


D 


1.9 x 10 41 cm" 2 


6.2 x 10 40 cm" 2 


6.2 x 10 40 cm" 2 


88% 


3 He 


1.9 x 10 41 cm" 2 


6.2 x 10 40 cm" 2 


6.2 x 10 40 cm" 2 


87.5% 



TABLE I. Integrated luminosities, and the effective polarization of the proton (neutron for D and 
3 He) in the projections for different ion beams and EIC energy configurations. 



generator [136 L The PDF parametrization used in EICDIS is accor ding to GRSV2000 
NLO model 1371 ]. while the input of PYTHIA is the standard scenario 1361 ] . We compare 
the SIDIS process for charged pion production between our results and those from the two 
models mentioned above for the 11+60 GeV EIC configuration in Fig. O and Fi g. ITU1 In this 
comparison, the fragmentation function by de Florian, Sassot, and Stratmann 13(j| is used. 



While the shape in the cross sections we obtain is in good agreement with those from the 
two models in the low Pt region, the results from the two models have been scaled down by 
a factor of 1.5 shown in Fig. |HJ Our results are about 1.5 lower than those from EICDIS and 
PYTHIA. Such a difference is likely due to missing contributions of longitudinal polarized 
photon as well as diffractive processes in our approach. Therefore, we take the conservative 
approach of using our rates for projections. On the other hand as shown in Fig. [TOj there 
are considerable differences between our results and those from the two models in the high 
Pt region. However, the differences become smaller in the higher Q 2 region (Q 2 >10 GeV 2 ). 
While more studies are needed in order to understand these differences, we decided to use 
correction factors to adjust the distributions from our approach to match those from the 
two models in order to make projections in the high Pt region (Fig. 15). At this stage of 
the study, this temporary solution is probably adequate. An EIC machine will naturally 
address this Pt dependence and provide the four- dimensional description of the unpolarized 
cross sections. 



3. Projections 

In this section, we present the projected results of TSSA at an EIC. Table U summarizes 
the used run time distribution, luminosities, and effective polarizations for different ion 
beams and energy configurations. In addition, we assume polarizations of ion beams to be 
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FIG. 9. The SIDIS comparisons at Pt <1 GeV/c for charged pion electroproduction. The EIC 
configuration is 10+60 GeV. 




Q 2 >10GeV 2 



This Work 
Pythia_standard 
EICDIS GRSV 2000 NLO 



6 7 
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FIG. 10. The SIDIS comparisons at Q 2 >1 GeV 2 and Q 2 >10 GeV 2 . 



70% and an overall detecting efficiency of 50%. The simulated data are binned according to 
different statistical precision for the TSSA measurement in different Q 2 regions. Different 
precisions are also chosen for different configurations. In particular, for both the 11+60 GeV 
and 11+100 GeV configurations, the statistical precision for each kinematic bin is set to be 
about 2.0 x 10~ 3 for Q 2 < 10 GeV 2 , and 4.0 x 10~ 3 for Q 2 > 10 GeV 2 . For the 3+20 GeV 
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FIG. 11. Collins/pretzelosity asymmetry projection with proton on ir + in a particular and z 
bin along with the calculated of asymmetries. The position of the dots are according to the Q 2 
axis on the left and the x axis, while the error bar of each dot is according to the scale of the 
asymmetry axis on the right. The calculated asymmetries are also according to the asymmetry 
axis. The black, green, and red dots represent the 11+60 GeV, 11+100 GeV, and 3+20 GeV EIC 
configuration. 



configuration, 4.0 x 10 3 , and 5.0 x 10 3 are chosen for Q 2 < 10 GeV 2 and Q 2 > 10 GeV 2 , 
respectively. 

Fig. [H] (Fig. [12]) shows the expected projection of 7r + Collins/pretzelosity (Sivers) asym- 
metry with a proton beam at a high luminosity EIC in the kinematic bin of 0.4 < z < 0.45 
and 0.4 GeV < Pt < 0.6 GeV. The x-axis represents Bjorken x, and the left y-axis is Q 2 . 
The position of each point in the plot represents the position of the kinematic bin in the 
x-Q 2 phase space. The right y-axis is the asymmetry. The error bar of each point follows the 



right axis. Together with the projection, several asymmetry calculations are a 
The codes to calculate the Collins and pretzelosi ty asymmetries are from 



107 



so p resented. 



Sivers asymmetry calculation is from 



139l | . and the 



571 ] . and 140| (red line). In the calculation, the PDF 
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FIG. 12. Sivers asymmetry projection with proton on ir + in a particular Py and z bin along with 
the calculated of asymmetries. The rest of the caption is the same as those in Fig. [TTJ 



is from MRST2004 parametrization 1411 ]. and the FF is from Ref. 1421 ] . Ref. 1071 ] pro vide s 
the Collins and Pretzelosity distributions, in which the Pt dependence is from Ref. 143 ] . 



The Sivers TMD is according to Ref. 



571 ] and the recent result of Anselmino et ai, and 



the Collins FF is according to Ref. 1431 ] . The calculated asymmetries also follow the right 
y-sixis of the plot. 

The selected 4-D projections for the average of the Sivers, Collins and pretzelosity asym- 
metries for the entire phase space are shown in Fig.[13]for tt + with a proton beam. The entire 
z coverage from 0.3-0.7 is divided into 8 bins (four z bins are shown in Fig. [T3|) ). We limit 
the projection at low region (Pj- < 1 GeV/c), where the Pt coverage from to 1 GeV/c 
is divided into 5 bins (three Pt bins are shown in Fig. [TBI) . In Fig. [T31 the central value of 
z bins increases from the left to the right. The central value of Pt bins increases from the 
top to the bottom. In addition to the proton results, the neutron results can be obtained 
with polarized 3 He and D beam. The selected 4-D projections of the corresponding neutron 
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FIG. 13. 4-D projections with proton on n + (0.3 < z < 0.7, GeV/c< Pt <1 GeV/c). The black, green, and red dots represent the 11+60 
GeV, 11+100 GeV, and 3+20 GeV EIC configuration. The rest of the caption is the same as Fig. [11] and Fig. [TZl 
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FIG. 14. 4-D projection with 3 He on ir + (0.3 < z < 0.7, GeV/c< < 1 GeV/c). The blue dots represent the projection from the 11-GeV 
SoLID SIDIS experiment, which overlap with the red dots (3+20 GeV). The rest of the caption is the same as Fig. [TT1 and Fig. [T2J 



results on 7r + using a polarized 3 He beam are sho wn i n Fig. [TUlj. Together, the projected 



ig. nsfl 



results for the 11 GeV SoLID SIDIS experiment |lll| are shown as blue points. The low 
energy EIC configuration (3+20 GeV) would provide the data connecting the phase space 
from both the fixed target experiment at 12 GeV JLab and the high energy 11+60 GeV EIC 
configuration. Furthermore, with additional kaon particle identification, the kaon SIDIS 
results can provide additional handle for the flavor separation, since kaon results would also 
tag the strange quark contribution from the sea. Fig. [T5l shows the projected results of K + 
on proton in the selected 4-D phase space. Since the kaon rates are normally about one 
order of magnitude lower than those of pions. The total number of points is significantly 
reduced. 

In addition, the high center-of-mass energy s at EIC would enable the studies of TSSA 
in high Pt region, where the twist-3 contribution will be large, and the intermediate Pt 
region, where one expects both the TMD and twist-3 formalism to work. In addition, a 
large Pt coverage of TSSA would provide the chance of forming Pt weighted asymmetry 
which is free of the Gaussian assumption of the transverse momentum dependence for both 
the TMDs and FFs. Fig. H~6l shows, as an example, the Pt dependence of the 4-D projection 
with proton on 7r + in one z bin. The number of points is limited at high Pt, where the 
differential cross section decreases. 

From these projections, it is clear that the TSSA can be precisely mapped in the full x, 
Q 2 , z and Pt 4-D phase space with a high luminosity EIC - a complete experiment (Table I) 
with a luminosity in excess of 10 34 cm~ 2 • will need 600 days of data taking. In particular, 
the EIC would facilitate the exploration of high Q 2 high x, and low Q 2 low x phase space. 
Furthermore, the large coverage of Pt would explore the TSSA in the high Pt region for 
the first time with SIDIS. The high luminosity is essential to realize the multi-dimensional 
mapping and extend the TSSA measurements to the unexplored regions (high Pt, high Q 2 
etc.). 

4- Projections for Pt- dependence of spin-azimuthal asymmetries 

Significantly higher, compared to JLabl2, Pt range accessible at EIC would allow for 
studies of transverse momentum dependence of different distribution and fragmentation 



There is a similar plot for Deuteron, which is not shown here. 
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FIG. 15. 4-D projection with proton on K + (0.3 < z < 0.7, GeV/c< Pt <1 GeV/c). The rest of the caption is the same as Fig. [TTJ 
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FIG. 16. 4-D projection with proton on ir + in one z bin to show the Pr dependence (0.45 < z < 0.50) in a range of 1 GeV/c< <10 
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FIG. 17. Double-spin asymmetry, A^, for positive pion production, using 4 GeV electrons and 60 
GeV protons (100 days at 10 34 cm 2 sec 1 ), as a func tion of Pt compared to published data from 



CLAS 



1451 ] (left). The right plot shows projections for 



471 ] and projected CLAS 12 measurements 
the same runnin g co nditions for the higher-twist lepton spin a sym metry compared to published 
data from CLAS [79] and HERMES Q and projected CLAS12 144] in one x, z bin (0.2 < x < 0.3, 
0.5 < z < 0.55). 



functions as well as transition from TMD regime to perturbative regime. Measurements 
of double-spin asymmetries as a function of the final hadron tra nsverse momentum at EIC 



will extend (see Fig. ITT]) measurements at JLAB12 144j, Il45| to significantly high er Pj 
and Q 2 allowing comparison with calculations performed in the perturbative limit 146 ]. 
Extending measurements of Pr-dependent observables to significantly lower x will provide 
access to transverse momentum dependence of quarks beyond the valence region. Much 
higher Q 2 range accessible at EIC would allow for studies of Q 2 -dependence of different higher 
twist SSAs, which, apart from providing important information on quark-gluon correlations 
are needed for understanding of possible corrections from higher twists to leading twist 
observables. 
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D. Transverse Single Spin Asymmetry for D (D) Mesons 



As discussed in Sec. I, the transverse single spin asymmetry of semi-inclusive neutral 
D meson production at large transverse momenta Pt would open a new window to the 
unexplored tri-gluon correlation functions. In this se ction , we discuss the potential of an 
EIC in this measurement using PYTHIA simulation 1471 ] . which has been compared with 
PEPSI in the inclusive electron scattering, and with PEPSI and our newly developed cross 



section weighting approach for the SIDIS pion production in Sec. C. 
projections with calculations based on an earlier paper on this subject 



iere, we compare our 



921 ] . New calculations 



based on the latest development 93j is ongoing. The mass of the charm quark is taken as 
1.65 GeV and the rest of the input parameters are from Ref. 1361 ] . The main channels of 
interest are 



D°(cu) ix + (ud)K- (su) 
D°(cu) ir-(ud)K + (us) 



(13) 
(14) 



where the branching ratio is 3.8 ± 0.1 % [148|. The two other main decay channels: D— > 
K + 7t~tt° and D— > K + 7t~7t~tt + are under investigation. 

The simulated phase space is limited in 0.05 < y < 0.9 and Q 2 > 1.0 GeV 2 for the 
11+60 EIC configuration with a proton beam. The main channel is the "direct" production 
channel, which is one of the four reaction mechanisms for producing neutral D mesons 
(vector meso n do minance (VMD), anomalous, direct and DIS) as shown in Fig. [TBI modeled 



in PYTHIA 



1491 ] for the hard 7 — p interactions. The Pt distributions of simulated D and 



D events are shown in Fig. dH] for all four mechanisms. The "direct" process dominates the 
production at Pt > 1 GeV, and the largest contamination is from the "anomalous" process. 
In the following studies, additional cuts of z > 0.15 and Pt > 1 GeV are applied. Fig. 120 
shows the momentum-polar angle distribution of electron and D meson in the lab frame. 
Fig. |2~T~1 shows the momentum-polar angle distribution of the 7r and K from the D (D) meson 
decay. The minimum momentum cut-off is chosen to be 0.6 GeV, and the minimum polar 
angle is chosen to be 10 degrees. 

In order to make projections, the D (D) mesons have to be reconstructed from the 
measured ir and K for each of the generated events. The signal-to-background ratio would 
strongly depend on the detector resolution. In this study, the momentum resolution is 
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a) VMD 



b) Direct 





c) anomalous 



d) DIS 



FIG. 18. Four diagrams of hard 7/? interactions modeled in PYTHIA. 



1 2 
log io (P T ) (GeV) 




1 2 

log io (P T ) (GeV) 



FIG. 19. Contributions from four interaction processes modeled in PYTHIA vs. Pt for D 
(l))-meson production. At large Pt, the "direct" production process dominates. The largest 
contamination is from the "anomalous" production at Pt > 1 GeV. 



assumed to be 0.8 ( , 



'0 • 



WGeV 



150j . The resolutions of polar and azimuthal angle are assumed 



to be 0.3 and 1 mr, respectively. The resulting resolution of the reconstructed invariant mass 
of the D meson is 1.8 MeV as shown in Fig. [221 In this case, the overall signal-to-background 
ratio is about 1.6 to 1. The background under the D meson invariant mass peak is due to 
the random coincidence of unrelated 7r and K in the final states after applying all cuts. The 
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FIG. 20. Top (bottom) two panels show the momentum vs. polar angle distribution for electron 
and D (D) in the lab frame. 

dilution of the background to the measured asymmetry is assumed to be: 



where S and B represent the signal and background within the invariant mass cut (1.86 < 
Mb < 1.87 GeV). The dilution due to the non-" direct" D meson production is applied in 
the same manner as in the final projection. 

Fig. [23] shows the projection for D (D) transverse SSA measurement for a running time 
of 144 days of proton beam at a luminosity of 3 x 10 34 /cm 2 ■ In this plot, the effect 
of 7r and K decay is included. The overall detection efficiency for this triple coincidence 
process is assumed to be 60%. The polarization of the proton beam is assumed to be 80%. 
An additional factor of \/2 Q is included in order to take into account the loss of precision 
in the angular separation. The data are binned 2-by-2 in terms of x and Q 2 . Within each 



7 Here, we adopt the first order approximation. As illustrated in previous section, such a factor, which is 
close to y/2, would depend on the azimuthal angular coverage and event distribution. 



33 



|D meson I |D meson 




P. (GeV) e„ (deg) 



FIG. 21. Left two panels show the momentum of tt vs. that of K from D (D) meson decay in the 
lab frame. Right two panels show the polar angle of tt vs. that of K. Red lines show the expected 
momentum and polar angle cut-off due to a realistic acceptance of an EIC detector. 




34 



0.1 



0.05 




0.1 



D meson S5A 144 days P L=3e34 on prolor 



if If If 



Q 2 ~1.3 GeV 2 
x- 0.0012 
P T ~1.6 GeV/c 



m 55A 144 days @ L=3e34 on protoi 



~| |i | | i 



if if if 



Q - 1 .3 GeV 
x ~ 0.0013 
P T ~1.6 GeV/c 



Q -1.3 GeV 
x - 0.0012 
z ~ 0.40 



D meson SSA144 days @ L^3e34 on proton B meson SSA 144 days @ L=3e34 on proton 



flUlJJJ, i L 

tf~1.3 GeV 2 — 



x~ 0.0013 
2 - 0.40 



0.05 




0.1 



P T -1.7 GeV/c 

~*^~\ 'j* 'l' 'I' ^ 'j 1 



Q 2 - 2.7 GeV 
x~ 0.0019 
P T -1.7 GeV/c 



Q - 2.8 GeV 2 
x - 0.0019 



Q 2 - 2.7 GeV 2 
x~ 0.0019 



0.05 




0.1 



Of ~ 2.2 GeV 2 





Q - 2.2 GeV 
x - 0.0059 
z - 0.38 



Q - 2.1 GeV 
x - 0.0061 
2 - 0.40 



0.05 



- 6.6 GeV 
x - 0.0082 
P T - 2.1 GeV/c 



Q ~ 6.4 GeV 
x ~ 0.0086 
P T - 2.1 GeV/c 




Q ~ 6.4 GeV 
x ~ 0.0086 
z - 0.37 




If If 



0.2 0.4 0.6 0.I 0.2 0.4 0.6 0.I 1 
z z 



4 1 
Pt 



Pt 



FIG. 23. The projected results on transverse SSA of D (D) meson production. The data are 
binned 2-by-2 in terms of x and Q 2 . Within each x-Q 2 bin, the projections are either plotted with 
z or Pt- The central kinematics are also listed in each plot. The first (second) column shows the 
projected results vs. z for D {D) meson. The third (fourth) column shows the projected results 
vs. Pt for D {D) meson. 



x-Q 2 bin, the projections are either plotted with z or Pt- The central kinematics are also 



listed in each of the panel. We show also theoretical predictions from Ref. 



92 



151]. 



E. Requirements on an EIC detector 

In this section, we summarize the requirements on an EIC detector from SIDIS processes. 

• Electron detection: 

As shown in Fig. [3], with the DIS cut of Q 2 > 1 GeV 2 , there is no need to cover the 
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extreme forward angle. The minimum polar angle coverage can be estimated by: 

min ~2axcsin(l/P e °), (16) 

where P e ° is the momentum of the incident electron beam. In addition, most of the 
scattered electrons have large momenta. Therefore, it is desirable to have: 



in order to achieve a good resolution of Bjorken x 150 ]. 



SIDIS 7T or K production at P T < 1 GeV: 

As shown in Fig. HJ leading hadrons span a large polar angular coverage in the lab 
frame. The momenta of most leading hadrons would be limited between 0.8 GeV 
and 7 GeV. Therefore, the separation of p/ir/K between 2.5 and 170 degrees and for 
momenta smaller than 7 GeV is essential to the success of the SIDIS program. In 
addition, a lower momentum cut-off for hadrons will enhance the overall acceptance. 

SIDIS 7T or K production at P T > 1 GeV: 

As shown in Fig. the high Pt events favor a large hadron momentum in the lab 
frame. Therefore, the separation of p/n/K for momenta larger than 7 GeV would be 
very useful for the high Pt SIDIS physics. 

SIDIS D or D production at P T > 1 GeV: 

As shown in Fig. [5H the momentum of the 7r and K (decay products of D meson) is 
in general smaller than 5 GeV. Therefore, a separation of p/n/K between 2.5 and 180 
degrees and for momenta smaller than 5 GeV is adequate for the identification of D or 
D meson for the transverse SSA physics. The more challenging requirement is on the 
detector resolution. As shown in Fig. [221 a 0.8% • 10 Q eV momentum resolution, a 0.3 
mrad polar and a 1 mrad azimuthal angular resolutions lead to a 1.8 MeV resolution 
of the reconstructed invariant mass of the D meson. A better detector resolution will 
lead to a better signal-to-noise ratio, which then leads to a better precision. 

Luminosity and Energy Coverage: 

Due to the mult i- dimensional nature of SIDIS processes, the 4-D (x, Q 2 , Pt and 
z) mapping of transverse SSA is essential for the success of TMD physics through 
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SIDIS at an EIC. In addition, from Fig. [5j it is essential to cover a few different CM. 
energies to reach ultimate mapping of asymmetries in x-Q 2 phase space. Therefore, it 
is essential to have L > 1 x 10 34 /cm 2 • s _1 at multiple s values (e.g. 11+100, 11+60, 
and 3+20 configurations). 

In this paper, we have summarized a recent workshop held at Duke University on Partonic 
Transverse Momentum in Hadrons: Quark Spin-Orbit Correlations and Quark-Gluon Inter- 
actions. The workshop participants identified the SSA measurements in SIDIS as a golden 
program to study TMDs in both the sea and valence quark regions as well as to study the 
role of gluons, with the Sivers asymmetry measurements as examples. A high-intensity EIC 
with a wide center-of-mass energy range will allow for studies of TMDs in multi-dimensional 
phase space with high precisions in both the valence quark region at high Q 2 and the unex- 
plored sea quark region, and allow for studies of tri-gluon correlation functions. Such studies 
will greatly advance our knowledge about the structure of the nucleon in three dimensions 
and transverse spin physics. 
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